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1. Introduction 

Recent years have seen significant developments in atom optics because of miniaturization, 
which has allowed workers to trap and manipulate atomic ensembles near a substrate via 
magnetic fields generated by currents in microfabricated wires jj-y]. These advances have 
permitted generation of Bose-Einstein condensates on atom chips, and production of coher- 
ent sources of atoms These scalable micro structures promise various applications on 
quantum information processing and micro-interferometry on a chip Q-Q. In both cases, 
an important step is to implement a local atom detector able to detect low atom numbers - 



ideally single atoms. Several detection possibilities have been, theoretical 
tally, demonstrated using the interaction of atoms with resonant light 11 



y and experimen- 



16| . For detection 



in an optical cavity, an important figure of merit is the "cooperativity" . The cooperativity 
C = <? 2 / K 7 describes the competition between the atom-field coupling g and the damping 
rates given by the half linewidth of the cavity and the atomic transitions (k, 7) respectively 



10Nl2|. In other words, C 1 is the number of atoms needed to significantly disturb the 



cavity. In a simple Fabry-Perot cavity, the cooperativity C oc J 7 / A depends mainly on the 
cavity parameters such as the finesse T and the intra-cavity beam cross section A. Even for 
a low cavity finesse, it is possible to reach useful values of the cooperativity if the intra-cavity 
cross section is small. In a system of several coupled cavities, as we will analyse here, the 
cooperativity is not a simple function of the finesse and beam cross section. This is because 
the electric field in different parts of the cavity can be very different and can depend on 
other parameters such as the lengths of the various sub-cavities. We therefore undertake in 
this paper to analyse a specific coupled cavity system with a view towards optimizing the 
cooperativity. We will show that the optimum operating conditions are not always easy to 
guess. 

In a prev.ous paper Q, we de 8 e ri bed a eav.ty consisting of an AlGaAs waveguide with 
mirrors at both ends. A cross section of the waveguide is shown in Fig. [1] (a). To use 
this device as an atom detector, we propose to open a small gap in the waveguide. Atoms 
within the gap can be detected or manipulated via their interaction with the intracavity 
field. Because of the considerable effective refractive index of AlGaAs (n ~ 3.1), the gap- 
vacuum interface has a reflectivity of 27 % and the gap thus forms a low finesse cavity 
within the larger waveguide cavity. The various length scales of the entire system are shown 
in Fig. [1] (b). In order to avoid excessive atom-surface interactions, it is necessary that the 
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gap width d be at least 2 fim 17j. On the other hand, the minimization of diffraction losses 
within the gap dictates that it be as narrow as possible. We have therefore concentrated 
on analyzing the characteristics of this system for gap sizes of about 2 /im. In our earlier 
work we gave an estimate of the cooperativity in the regime where the finesse is limited 
by the waveguide propagation losses. In the present paper we discuss a more systematic 
investigation. First we study numerically and analytically the transmission and reflection of 
the gap as a function of its width, taking into account the losses due to diffraction. Knowing 
the transfer matrix of the gap, we estimate in section [3] the finesse of a resonator based on 
such a opened waveguide. In section HI we compute the amplitude of the field in each part 
of the resonator in order to calculate the coupling parameter g of an atom to the light field. 
These calculations allow us to calculate the cooperativity and optimize it as a function of 
the cavity parameters. 

We confirm the approximate results of Ref. 17|. However, we believe it should be possible 
to fabricate a waveguide with sufficiently low loss that the total losses are dominated by the 
unavoidable loss in the gap. In this case, we show that the optimization of the cavity is very 
different. 

2. Transmission and reflection of the gap 

We begin by giving some more details of the calculations of the losses induced by the gap 
which allowed us to estimate the cooperativity in our previous paper 17|. To perform 
our analysis, a key issue is the spatial variation of the eigenmode of the waveguide as it 
propagates in free-space. Because of the complex geometry of the waveguide (Fig. [TJ (a)), 
the spatial mode profile must be calculated numerically. We consider a field that propagates 
inside the waveguide for z < 0, and that reaches at z = a waveguide/vacuum interface 
(Fig. [U (b)). Using a 3 dimensional fully- vectorial aperiodic- Fourier modal method 181 ] . we 
compute the complex amplitude £ (x,y, z) of the field that propagates in free-space in the 
half-space corresponding to z > 0. Numerically, we perform the simulation from z = to 
z = 35 /im. We define the "axis of the field" as the coordinates (xo;2/o) that correspond 
to the maximum of \£o(x, y, 0)| and set the origin of the coordinates so that x = Uo — 0- 
We define E = £o(0, 0, 0), and then rewrite £o(x, y, z) = E f(x, y, z), where / is the mode 
profile during it's propagation along the z-axis. Due to the continuity of the electric field, 
f(x, y, 0) corresponds to the profile of the eigenmode inside the waveguide. We can therefore 
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compute the overlap between the propagated field and the field inside the waveguide, which 
is expressed as: 

Q{z) = (1) 



/ \f(x,y,0)\ 2 dxdy J \f(x,y,z)\ 2 dxdy 
For z > 0, the overlap function is |Q(z)| < 1, which expresses the fact that because of the 
diffraction during the propagation in free-space, the field cannot be totally recoupled into 
the waveguide on the other side of the gap. The knowledge of f(x,y,z) also allows us to 
compute the phase and amplitude of the electric field on the axis of the waveguide as the 
field propagates: 

/(0, 0, z) = |/(0, 0, z) | exp i(kz + <p (z)), (2) 

where ifo is the correction to the plane wave phase and varies slowly with zona scale A. 
Once we know the free-space propagation of the eigenmode, we can calculate the transmission 
and reflection of a gap of width d. We consider the gap as two plane interfaces separated 
by a distance d. Each interface has a transmission and a reflection given by the Fresnel 
coefficient of reflection and transmission between a material of refractive index n — 3.15 and 
the vacuum no = 1. On the left interface, the light coupled into the gap has the spatial 
profile of the eigenmode of the waveguide. The mode inside the cavity is considered as the 
superposition of £+ and £~, the fields corresponding to the 2p th (propagating from left to 
right) and 2p + 1 th (propagating from right to left) reflections at the gap interfaces. The 
amplitude {S~) can be expressed as a function of the incident field amplitude E t as 
Ep = tr 2p Ei and E~ = tr 2p+1 Ei. The Fresnel coefficients t and r correspond respectively to 
the transmission coefficient for a wave travelling from the waveguide into the vacuum, and 
the reflection coefficient of the vacuum/ waveguide interface. The spatial profile of and 
£~ is related to f(x, y, z) by: 

£p(x,y,z) = E+f(x,y,2pd + z) (3) 
£;{x,y,z) = E;f(x,y,{2p + 2)d-z). (4) 

A fraction t' , where t' is the transmission coefficient for a wave travelling from the gap 
into the waveguide, of each component is transmitted through the right interface, but 
since the light diffracts during its propagation in the free-space of the gap, only a fraction 
Qp — Q(2pd + d) of the transmitted field is coupled back into the eigenmode of the right 
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waveguide. The amplitude E t of field transmitted through the gap is therefore: 

E t = t'J2Qptr 2p E t . (5) 

p 

In a similar way, one can write the amplitude of the field reflected by the gap as the sum of 
the amplitude of the light directly reflected by the waveguide/vacuum interface — rEi and 
the mode-matched part of £~ that is transmitted back into the left part of the waveguide: 

E r = -rEi + t'J2 Q;tr 2p+1 E U (6) 
p 

where Q~ = Q(2pd + 2d). 

Considering the gap as a single entity, its complex reflection and transmission coefficients 
for the amplitude are r g = E r /Ei and t g = E t /Ei, and characterize completely the gap. 
The phase and amplitude of f g and t g are shown in Fig. [2j To confirm our results, we can 
also use the simulation software to directly simulate the propagation through a waveguide 
- vacuum - waveguide stack. We denote by r g ^ sim and t g)Sim the reflection and transmission 
coefficients for the amplitude deduced from the numerical simulations. Figures [2] (a) and [2] 
(b) compare the result of the analytical model to the numerical simulations. We see that the 
amplitude and phase match very well. We also compared the losses, defined as 1 — |t s | 2 — |r 5 | 2 , 
calculated using the two methods and they agree within 20 %. 

3. Resonator losses induced by the gap 

The losses calculated above only correspond to the losses for a single pass inside the cavity. 
The losses induced by the gap in a cavity will depend on the interferences between the two 
propagating fields that create the standing wave inside the resonator. Therefore, to calculate 
the maximum finesse of a cavity based on a gapped waveguide, we need to estimate the losses 
induced by the gap in the presence of the counter-propagating field reflected by a perfect 
mirror that terminates the cavity (reflectivity r = —1) (Fig. [T] (b)). We denote by E\ 
(E3) the amplitude just before (after) the gap of the field propagating in the left (right) 
waveguide from the left to the right, and E-x (E-3) the amplitude just before (after) the 
gap of the field in the left (right) waveguide propagating from right to left. Following the 
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notation of section [21 one can write: 



E-3 



E 



E 3 



t g {d)E 1 +r g {d)E_ 3 
t g (d)E^ 3 + f g (d)E 1 
r exp(2inkL 3 )E 3 , 



(?) 
(8) 
(9) 



where n is the effective refractive index of the eigenmode of the waveguide and k = 2tt/X 
the wavenumber of the light in vacuum. Solving (J7HH]), we can express E-i as a function of 



where r e ff is the effective reflectivity of the ensemble gap + waveguide + end mirror. Once 
we have calculated r e ff, we can calculate the finesse of a resonator with a gap. It is the 
finesse of the equivalent resonator made of a perfect mirror (reflectivity tq = —1) facing a 
mirror of reflectivity r e ff, separated by a distance L\. If we write f e // = r e ff exp(2z</>), we 
can express the finesse as: 



In the following, we assume that L 3 3> d. Since we are looking for a resonance of the 
cavity, we consider the response r e ff for wavenumber k + 5k such that 5k < l/L 3 <C 1/d. 
Therefore we will consider that r g (k + 5k) rs f g (k) and t g (k + 5k) ps t g (k) and that the only 
/c-dependence of <p comes from the exp(2inkL 3 ) of Eq. [9j Knowing the effective reflectivity 
f e ff, we are able to compute the finesse of such a compound micro cavity. The results are 
shown in Fig. [3] for different values of d, L\ and L 3 . We varied d between 1.5 and 2.5 /xm, 
and L 3 between lOOOA/n and lOOlA/n. Once d and L 3 are fixed, the distance L\ is no longer 
a free parameter since it has to be chosen to verify the resonance condition of the cavity. 
However, we can always choose L\ within X/2n of 1000A/n. Therefore, in the calculations, 
we assume that L\ lOOOA/n, and only plot the finesse as function of d and L 3 . The 
choice of lOOOA/n is arbitrary and simply ensures the two conditions L 3 ^> d and L\ ^> X/n, 
necessary to neglect the dependence of L\ with L 3 and d in the calculations. Since L 3 only 
appears in equations (JZHSJ) as exp(2iknL 3 ), one can check that 



E-i = r e f f (d,k,L 3 )E : 



(10) 




(11) 




(12) 
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where 5L% = L3 mod X/2n. Equation [12] clearly shows that J- a does not change if we 
multiply Li and L 3 by the same factor, as long as we keep SL 3 constant. 



The plots in Fig. [3] show that the finesse globally decreases when the gap width d 
increases. This is consistent with the result of Fig. [21 and corresponds again to the fact 
that the diffraction losses increase with d. In addition to the global decrease, we observe 
oscillations of the finesse as a function of both d and L3. To understand these oscillations, 
we recall that the losses correspond to fields that are transmitted across the gap but are 
not coupled into the guided mode on the other side. In the presence of a perfect right-end 
mirror, the amplitude of the field leaking out at each interface is the sum of a contribution 
proportional to E\ (the single pass losses of the incident field) and a contribution proportional 
to i?_3 (the single pass losses of the field reflected by the end mirror). Depending on the 
relative phase of E\ and E_ 3 , the interference between these two contributions can be either 
constructive or destructive, and one can reduce the losses in the gap and reach a much larger 
finesse than what could be deduced from the single pass losses. 

4. Atom-field coupling 

To estimate the cooperativity of the system, we need to calculate the amplitude of a one- 
photon field in the center of the gap. Following the same argument as in the section (3j the 
field in the gap E g (z) is the sum of the multiple reflections of the light coming from the left 
waveguide and the multiple reflections coming from the right. Consequently, using equation 
[21 one can rewrite E g (z) as the sum of E£(z) and E~(z) where 



E+(z) = [Ei tr p \f (pd + z)\exp(ikpd + i<p (pd + z)) 



p even 



+i?_3 tr p \fo(pd + z) \ exp(ikpd + itpo(pd + z))] exp(ikz) 



p odd 



= Eg{z) exp(ikz) 
E~(z) = [Ei tr p \f (pd + d — z)\ exp(ikpd + i(fo(pd + d — z)) 



(13) 



p odd 



+£L 3 tr p \f (pd + d — z)\ exp(ikpd + i(fo(pd + d — z))] exp(— ikz) 



p even 



= E g (z) exp(-ikz), 



(14) 
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where |/o(z)| = 1/(0, 0, z)\ is the amplitude of the electric field on the axis of the waveguide. 
Here, E£ and E~ vary slowly on the scale of A. As a result, the field inside the gap is the 
sum of two terms evolving at exp(ikz) and exp(— ikz). The interference of these two terms 
will not necessarily be constructive at the exact center of the gap z — d/2, but there will be 
an anti-node at a position z* within A/2 of the center where the two amplitudes will add up 
constructively and 



\E g (z*)\ = \E+(z*)\ + \E;(z*)\ « \E+ ( i\ | + \E~ (i ) |. (15) 



Up to this point, the amplitude of the field in the cavity has been set by an arbitrary 
factor Ei. To calculate the single-photon amplitude, one has to write that the total energy 
stored in the cavity is equal to Hcu. Since d <C {Li, L 3 }, we can neglect the energy stored in 
the gap. The total energy of the cavity is therefore 

hu ^ An 2 e A(Li\Ei\ 2 + L 3 \E 3 \ 2 ), (16) 

where v = c/A is the photon frequency, eo is the electric permittivity of vacuum, and we have 
assumed that \E-\\ \Ei\ since the losses are small. Finally, the single photon amplitude 
is given by: 

E = V^(\Et\ + I^D (17) 
^/An\ A(Li\Ei\* + L 3 \Etf) 

Assuming that Li m L 3 m L, we can rewrite the maximum atom field coupling g at the 
center anti-node as g = gog(d, L\, L 3 ) with: 



90 = \l TT-TTT, (18) 



k 2 n 2 AL 
and 

(l9) 

Here again, Li is chosen to ensure that the cavity is resonant at the atomic frequency, and 
is no longer a free parameter once d and L 3 are fixed. We can therefore plot g(d, Li, L 3 ) = 
g(d, L 3 ) as function of d and L 3 (Fig. H]). The value of g varies between 1 and n. This can 
be very well explained by the inset of the Fig. H] showing the electric field distribution inside 
the cavity. If Li, L 3 and d all correspond to an integer number of half- wavelength of the 
field, the interferences of the partial reflections at the air-waveguide interfaces enhance the 
field in the gap by a factor n, leading to g — n. If d corresponds to a half integer number 
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of half wavelengths in the gap, g = \j2n / \Jn 2 + 1 is constant. If d is an integer number 
of half wavelengths, but L\ and L3 are equal to A/4n (mod X/2n), then the amplitude of 
the standing wave of the electromagnetic field is constant in the whole cavity and g = 1. 
To get a feeling for the magnitude of the coupling, we use the parameters of the above 
cavity including the effective mode area A = 9.9 fim 2 (see Ref. The formula ( fl8]) gives 

go w 2ir x 30 MHz and so the maximum value of the coupling is about 2ir x 100 MHz. 

5. Cooperativity of the atom-cavity system 

The figure of merit of the cavity is given by the cooperativity C = g 2 / 'kj. The cavity loss 
rate k is inversely proportional to the finesse: 

K = 2n(L 1 + L 3 )T> (20) 



and the cooperativity is then given by: 



C = 2^-~g 2 (d,L 3 ), (21) 
TmA 

where a a b s = Qir/k 2 is the resonant atomic cross section. 

We can see that the cooperativity is proportional to the product g 2 ^. Comparing the 
results of sections II and III shows that optimizing the cooperativity will be a trade off 
between maximizing J 7 and maximizing g. The total finesse of the cavity is given by J 7-1 = 
J-'J 1 +J C ~ 1 , where corresponds to the finesse induced by the diffraction losses (calculated 
in section 3), and T v corresponds to the other intrinsic losses of the cavity (such as the losses 
due to the absorption or scattering of the light during its propagation in the waveguide part 
of the cavity). Using our numerical model, we have computed the cooperativity as a function 
of d and L 3 in the two regime 3> T p and Jrf < J p . Each time we vary d between 1.5 and 
2.5 jum, and L 3 between lOOOA/n and lOOlA/n keeping Li w 1000A/n in the two different 
cases. 

For the first case, we consider the losses due to absorption and/ or diffusion in the AlGaAs 
waveguide as experimentally measured in [17|], corresponding to T v ~ 92. The results are 
shown Fig. [5] and they confirm what was suggested in • If the propagation losses are 
high, the total finesse is limited by T p , and the maximum of the cooperativity is a point in 
the parameter space very close to the values d = 1.95 fim and nL 3 = 1000A which optimize 
the coupling g. Still, the calculation shows that the cooperativity can be increased by ~10 % 
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by slightly decreasing d and increasing L3, which slightly increases the finesse without losing 
too much of the coupling. 

We also performed simulations for a lossless waveguide, that is, in the situation where the 
only loss is diffraction in the gap (J-^ <C J-" p ). In this case the situation is completely different 
(see Fig. E]). Optimizing the cooperativity corresponds to maximizing the finesse rather than 
to maximizing the coupling. To illustrate this point, consider a gap width d = 1.95 /im. For 
77X3 = 1000. 25A, g 2 is n 2 ~10 times smaller than for 77X3 = 1000A. However, the finesse is 
~29 times higher for 77X3 = 1000. 25A and this choice clearly results in a higher C. By setting 
d = 2.05 /im and nL 3 = 1000. 23A, one can further increase the finesse and the coupling, 
reaching a cooperativity C max ~ 9.4. 

It is also interesting to notice that the maximum of C is very narrow with respect to L 3 , 
but is much wider with respect to d. In the vicinity of d ~ 1.95 /im, the C > 9 area is 180 
nm wide along the d axis. This is important experimentally, because unlike the waveguide 
parts of the cavity (L\ and L3), whose refractive index can be changed by changing the 
temperature or by injecting a current, the gap's optical length cannot be fine-tuned once 
the cavity has been fabricated. These calculations show that if the gap is fabricated within 
90 nm of the optimal width, it will always be possible to reach C > 95 %C max by slightly 
tuning the index of refraction of L3. 

Another important parameter characterizing the atom-cavity system is the saturation 
photon number jjj], the photon occupation number in the cavity necessary to saturate the 
atomic transition, jV sa t = (7/ ' g) 2 1Q|. The calculations in Sec.|4]show that typically the atom 
field coupling is at least 10 times the natural linewidth of the Rb transition (7 = 2-k x 3 MHz). 
Thus the saturation photon number can vary between 0.01 and 0.001 as the gap width and 
the total cavity length are varied. 

The saturation photon number is useful to characterize the cavity transmission in the 
presence of an atom. If we adjust the cavity parameters to give maximum cooperativity, the 
coupling is minimal (g ~ 2n x 30 MHz) while the finesse is about 1500. In this situation we 
have n ~ In x 34 MHz close to the value of g, but much larger than the atomic damping rate. 
Without an atom the cavity, the transmission spectrum close to resonance show a peak with 
a width of order k. Over a broader spectrum of frequencies, the cavity exhibits a complex 
structure characteristic of a multimirror Fabry-Perot interferometer as analyzed in Ref. 19 ]. 



The insertion of an atom into the cavity causes a dip in the in the absorption spectrum 
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when the light frequency is close to the atomic resonance [121 ] . The width of this dip is of 



order 7. Reference 



12] discusses the optimization of the signal to noise in such a situation 



and points out that the optimum signal to noise ratio is achieved when the photon number 
is equal the saturation photon number. Direct observation of the vacuum Rabi splitting is 
not possible because of the large cavity damping. 

6. Summary 

We have analyzed the cooperativity of a cavity consisting of an optical waveguide in which 
a gap is opened to let atoms interact with the cavity field. Because of the high index of 
refraction of the guided parts of the cavity, the Fresnel reflections at the interfaces are not 
negligible and the gap behaves like a low finesse Fabry-Perot inside the high finesse cavity. 
We have separately calculated the losses induced by the diffraction of the light in the gap, 
and the value of the atom-light coupling. Choosing the cavity parameters such that all sub- 
cavities are resonant with the electromagnetic field maximizes the atom field coupling, but 
this choice also maximizes the diffraction losses introduced by the gap. This choice is optimal 
as long as the losses are dominated by other losses, such as absorption or diffusion in the 
waveguide. If the non-diffractive losses are small enough, a higher cooperativity is obtained 
by choosing the gap size and cavity length to decrease the coupling strength and increase 
the finesse. Our simulation predicts that a finesse on the order of 1500 and cooperativity 
C ~ 9 should be possible. We also showed that the value of the cavity parameters, especially 
the gap width, are tolerant to small deviations. The possibility of tuning the waveguide 
lengths means that close to optimum performance can be obtained despite small errors in 
the gap width. The achievable cooperativity can attain values which should be of interest 
for the detection of single atoms on chips. Although the fabrication of a gap remains a great 
technical challenge and has not yet been achieved, the fully integrated nature of the device 
we have analyzed means that we can envision a system of several such structures fabricated 
in parallel on the same substrate. This would permit a structure similar to that in reference 



161 ] with the additional feature of cavity enhanced atom sensitivity. 
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Fig. 1. (a) Transverse geometry of the waveguide. The waveguide is a rectangular ridge 
of AlGaAs, 4 /im wide. The false color image shows the result of a numerical calculation 
of the amplitude of fundamental mode of the waveguide operating at A =780 nm. The 
mode area is calculated numerically and we find A = 9.9 /im 2 . (b) Top view of the gapped 
waveguide. Atoms are to be guided into this gap to interact with the intracavity field. The 
gap constitutes a low finesse cavity within a larger waveguide cavity which is closed by 
reflecting surfaces at the both ends of the waveguide. 
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Fig. 2. (a) Amplitude of the transmission and reflection of the g function of the gap 

width, (b) Phase of the transmission and reflection coefficients (the propagating term kod has 
been subtracted). The solid lines correspond to the analytical model (f g , t g ), the square and 
circles correspond to the result of the full numerical simulation (r g ^ im and t ffiS j m respectively, 
see text), (c) Losses induced by the gap. The solid line correspond to 1 — |t 3 | 2 — |r 9 | 2 , the 
dotted line + squares correspond to 1 — \t gtSim \ 2 — \ f g ^ im \ 2 . The losses are maximal for a gap 
width equal to an integer number of half wavelengths. 
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Fig. 3. (a) Cavity finesse as function the gap width d and waveguide length L 3 . (b) Finesse 
as a function of L 3 for d — 1.95 fim (solid red line), and d ~ 2.15 /im (dotted black line), (c) 
Finesse as a function of d for nL 3 /\ = 1000.5 (solid red line), and nL 3 /X = 1000.25 (dotted 
black line). 
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Fig. 4. Contour plot of the atom- field coupling g as function d and L3. Insets show the 
electric field profile inside the coupled cavity for special cases of d and L3 with a normalisation 
I Ei I + \ E-\\ = 1. (a) Electric field profile for d = 1.95 /xm and L 3 equal to an integer number 
of half wavelengths, (b) Electric field profile for d = 1.95 fim and L3 = A/4n (mod \/2n). (c) 
Electric field profile for d = 1.75 fim and L3 equal to an integer number of half wavelengths, 
(d) Electric field profile for d = 2.05 /xm and L 3 = 0.23A (mod A/2n) (value that optimizes 
C in the lossless case). 17 
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Fig. 5. Cooperativity of a Rb atom in the coupled cavity system as a function of d and L3, 
including propagation losses in the waveguides. The attenuation coefficient corresponds to 

17j. The coupling is maximized for d = 1.95 fim and nL 3 = 1000. 5A. 



that measured in Ref. 
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